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Abstract 

We investigate the possibility for a quark-antiquark pair to form a bound state at temperatures 
higher than the critical one (T > T c ), thus after deconfinement. Our main goal is to find analytical 
criteria constraining the existence of such mesons. Our formalism relies on a Schrodinger equation 
for which we study the physical consequences of both using the free energy and the internal energy 
as potential term, assuming a widely accepted temperature-dependent Yukawa form for the free 
energy and a recently proposed nonperturbative form for the screening mass. We show that using 
the free energy only allows for the IS bottomonium to be bound above T c , with a dissociation 
temperature around 1.5 x T c . The situation is very different with the internal energy, where we 
show that no bound states at all can exist in the deconfined phase. But, in this last case, quasi- 
bound states could be present at higher temperatures because of a positive barrier appearing in 
the potential. 
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I. INTRODUCTION 



The anomalous suppression of heavy meson production in heavy ion collisions has been 
proposed as a possible signature of deconfmement for about twenty years The basic 
idea is that the heavy mesons produced before the formation of a thermalized quark-gluon- 
plasma would tend to dissociate in the deconfmed medium because of the screening of the 
quark- ant iquark (qq) interaction by the large number of color charges in the medium. The 
mechanism is clearly analogous to the Debye screening by electromagnetic charges in QED. 
Since the pioneering work of Ref. the suppression of heavy quarkonium production at 
finite temperature - that is T > 0, T = corresponding to "usual" QCD - has been 
intensively studied in the literature (see Ref. for a review). In particular, the dissociation 
temperature of a particular meson, that is the temperature at which the quark and the 
antiquark become unbound, is particularly relevant in order to understand the mechanism 
of quarkonium dissociation in a quark-gluon-plasma. 

It is a well-known fact that potential models are able to accurately reproduce the ex- 
perimental meson mass spectra at T = !|. Basically, a meson is then seen as a qq pair 
interacting via a so-called Cornell potential ar — Aa s /3r (or any other QCD- inspired po- 
tential), r being the distance between the quark and the antiquark. The Cornell potential, 
validated by lattice QCD computations of the energy of a static qq pair 4j, contains a con- 
fining linear part and a Coulomb-like term which comes from one-gluon-exchange diagrams. 
In this potential, a can be interpreted as the tension of a flux tube of length r linking the 
quark and the antiquark, and a s is the strong coupling constant. Because of the success of 
potential models at zero temperature and of the idea that color screening implies modifica- 
tions of the Cornell interaction, several attempts to understand meson properties at finite 
temperature within the framework of a Schrodinger equation with a phenomenologically 



modified potential have been made 
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The most accurate way to find a relevant potential term at T > is again provided 
by lattice QCD simulations, from which the free energy between a static qq pair at finite 

, n 

temperature can be computed [9|, [10|, LL3J • A general result of lattice QCD calculations is the 
observation of a critical temperature, denoted as T c , above which the confining part of the 
free energy vanishes: That is the deconfmement (see for example Ref. 12[). The free energy 
obtained in these calculations, depending both on T and on r, can be used as a starting 
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point to build a potential model at T ^ 0. However, in some other works, it is argued 



that the internal energy should better be used as potential term 14|. 



161 ]. because of 



the entropy contribution which is present in the free energy. The situation of how should 



potential models be applied to qq states at T > is thus still not completely clarified [16|. 

We propose in the present work to investigate the existence of mesons above the critical 
temperature by using a Schrodinger equation with an appropriate potential term. The 
motivation to study such a temperature domain does not only come from intrinsic theoretical 
interest, but also because the current temperatures reached by experiments are in the typical 
range (1 — 2) x T c 13]. Let us note that what we call a meson in the following is a qq bound 
state, that is a state with a negative binding energy and an infinite lifetime - at least 
formally. Our model, that we present in Sec. [Ill is a rather simple one, but it contains the 
main qualitative features of most of the potential models which were previously developed. 
Moreover, the Yukawa form that we use for the free energy is in agreement with recent 
and accurate lattice QCD results ll|. In this work, we are mainly interested in qualitative 
and mostly analytical results constraining the existence of mesons at finite temperature. 
Actually we propose a simple method which can also be applied to a more complete version 
of the model studied here. In particular, we focus on the dissociation temperature, which is 
a physically relevant observable, leading to a direct picture of the evolution of the number 
of bound states with the temperature. We show that the determination of the dissociation 
temperature is eventually ruled by a single dimensionless parameter which corresponds to 
the strength of the potential in a rescaled set of coordinates. Constrains on the dissociation 
temperature are thus obtained from constrains on the strength of the potential such as it 
is attractive enough to admit bound states. Since two choices are generally assumed in 
the literature for the potential term, we study them both in Sees. IIIII (free energy) and IIVI 
(internal energy). As we will see in the latter, both possibilities lead to rather different 
predictions; we then summarize our results in Sec. |V] 
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II. THE MODEL 



The binding energy of a qq pair in a medium at finite temperature T > can be obtained 
by solving the following Schrodinger equation 



f + V(r,T) 



i,(r,T)=E rd {T)i,{r,T), (1) 

where /i = m q mq/(m q + irtq) is the reduced mass of the system (we work in units where 
h = c = 1). As it is generally done, we assumed a central form for the temperature-dependent 
interaction potential V(r, T). Consequently, the binding energy e n i (T) does not only depend 
on the temperature, but also on the radial quantum number n and on the orbital angular 
momentum i. An important question is the following: Should the free energy of the qq pair 
be used as potential term, or should the internal energy be used? Up to our knowledge, 
the answer is far from being unanimously accepted. Consequently, we will consider both 
possibilities in this work, and study their consequences on the existence of bound states from 
the analysis of Eq. ([I]) . 

It is known for a long time that lattice calculations can directly compute the free energy, 
denoted as Fi(r,T), between a quark and an antiquark placed in a thermal bath of gluons 
and light quarks 9j. For T > T c , recent lattice computations agree with a Yukawa form for 
the free energy of a qq pair in a color singlet, which can be parametrized by 11] 



A p -m D (T)r 

F 1 (r,T) = --a s (T)——, (2) 

where a s (T) is an effective strong coupling constant depending on the temperature. Equa- 
tion ([2]) shows that the free energy is screened in the finite temperature medium, in analogy 
to what happens in QED. The potential at T > T c is then no longer a confining one. The 
screening parameter could be considered as a QCD Debye mass, denoted as rri£,(T). Since 
the free energy naturally emerges from lattice QCD computations, many works directly 
assumed that V(r,T) = F x (r,T) to study the properties of mesons with the eigenequa- 
tion ([!]) 0, [lsj. We point our that, if these works use the free energy as potential term, they 
do not necessarily take the same expression as ours for Fi(r, T). Actually, the simple form 



(T5]) has the advantage of being in agreement with the recent results of Ref. [11). 

On the other hand, the free energy contains an entropy contribution, and other approaches 
have suggested that the internal energy, defined as usual in thermodynamics as 

U 1 (r,T)=F 1 (r, T) + T S 1 (r,T) = F 1 {r,T) + T (r, T) , (3) 
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Ml 



should be preferentially used as potential term in the Schrodinger equation ([I]) {2). 
The reason for such a choice could be related to the time scales involved in the system 
The first time scale that should be considered is r& is the typical time associated to a 
particular bound state, i.e. (r/r). The second one is r^, that is the time which is needed to 
transfer heat to matter by changing the entropy Si. If T& <C r^, then the heat transfer can 
be neglected, and the internal energy is the relevant potential term. If not, the free energy 
should be used. With the free energy defined as ([2j), the internal energy reads 

-m D {T)r 



U t (r,T) = ~ [(Ta s (T))' - T a s {T)m D {T)'r] 6 



(4) 



where the prime denotes a derivation with respect to T. 

The evolution of Fx (r, T) and Ux (r, T) with the temperature is completely fixed by the 
functions a s (T) and mo{T). For the running of the strong couplin g c onstant with the 
temperature scale, we will assume the well-known one-loop expression 10] 

2tt 



a s (T) = 

(11 - |JV» In 

From lattice computations, we choose the value 

A„ = PT C , with = 0.104 ±0.009, 



(5) 



(6) 



which has been obtained in Ref. [lOj. Actually, a s {T) is also known at two-loops [lOj, but the 
one-loop formula already captures the essential physical features of the two-loops running 
coupling constant. Since we are mainly interested in a qualitative description of the existence 
of mesons versus the temperature, formula (jSJ) is thus sufficient for our purpose. Lattice 
QCD also provides a critical temperature T c appearing in a s (T). Following the number of 
light flavors Nf, different values can be found, which globally lie in the range 150 — 300 



MeV 08 



201 ] . We will here take the recent estimation of Ref. 19j, where Nf = 3 with 



two light quar ks of the same mass (u and d) and one heavier (s). It is computed in this 
work that [19] 

T c = 169 ± 16 MeV => A CT = 17.6 ± 3.2 MeV. (7) 



If the function giving a s (T) is generally acce pte d, severa 



in the literature for the screening mass m£>(T) 
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21 



22 



different forms can be found 
23j, depending mainly on the 



temperature range which is considered. In this paper, we will assume a recently proposed 



nonperturbative formula, obtained in Ref. 21] within the framework of the background 
perturbation theory, which states that 



m D (T) = 47r?7 c a a s (T) T, 



(8) 



with c CT = 0.566 ± 0.013 (lOj and rj = 2.06 [2l|. Defining 7 = lirqc, we have 



mn(T) = 7 a s {T) T, with 7 = 14.652 ± 0.337. 



(9) 



Contrarily to the well-known form m D (T) = [4tt (1 + Nj/6)] 1 / 2 [^(T)] 1 / 2 T which is 
assumed to hold at T ^> T c only [23j, formula ([9]) is expected to hold also at T > T c . As we 



show in Fig. [TJ formula fits indeed rather well the lattice QCD data, taken from Ref. 11], 



in the range (1 — 3) x T c . We also tried to fit the data with a perturbative-inspired ansatz 



of the form 



(T) = A f y/4ir{l + N f /6)yfe(X)T. 



(10) 



Even with the ad hoc value Af = 1.65, which leads to the best fit of the data, formula ([10]) 
clearly yields a poorer agreement with lattice QCD (see Fig. CD). Consequently, it is better 
justified to work with the screening mass defined by Eq. Q). 

Thanks to Eqs. (jSJ), (CJ) and ([9]), the free and internal energies are completely known, and 
it is possible to compute the binding energies e n £ (T) from the eigenequation (flj). Then, the 
dissociation temperature of a particular [(n + l),^o] state (in spectroscopic notation) can 
be computed, that is the temperature T no i for which 



(11) 



The physical meaning of such a temperature is the following: Once T > T noio , no qq bound 
state with quantum numbers n > or I > can exist. 



III. FREE ENERGY AS POTENTIAL TERM 



Let us now investigate the existence of bound states when the free energy ([2]) is used in 
the Schrodinger equation (OQ): V(r, T) = Fi(r, T). For further computations, it is convenient 
to work with the dimensionless variables x, q, defined as 

x = mr,{T) r, q = m£>(T) _1 p. (12) 
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FIG. 1: Graph of mr>(T)/T as a function of T/T r for a qq pair in a color singlet. The full circles 



are taken from the lattice computations of Ref. 



XII ] . Formula Q, computed with the mean values 



of the parameters, is plotted with a solid line. The error bars on the parameters appearing in 
mc(T) actually allow to span the shaded area. The perturbative-like formula (|10p is also plotted 
for comparison with the fitted value Af = 1.65 and the mean values of the parameters (dashed 
line). 



Equation (JT]) can then be rewritten as 



[q 2 - g P(x)] ip(x, g) = E nl (g) ip(x, g), 



with 



and 



P(x) 



x 



£m (g) = 

\/i a s (T) 



2/1 £ n£ (T) 

m D (Tf 



> 0. 



(13) 



(14) 



(15) 



3 m D (T) 

Equations (I13H15I) show that the dimensionless parameter g ruling the depth of the 
Yukawa potential, —gP(x), is the only physical parameter which is responsible for the 
existence of bound states. In particular, the numbers g n i such that £{g n e) = are the 
critical strengths of the potential: When g < g no e , no bound states with quantum numbers 
n > uq or i > £q are present in the spectrum. The numbers g n i actually only depend on the 
form the potential, which is —ge~ x /x in this section. They are given in Table [I] for some 
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TABLE I: Critical strengths g n £ of the potential — ge x /x for some values of n and I. 



n/£ 


1 





1.680 9.082 


1 


6.447 17.745 



values of n and £. The lowest value, namely goo = 1.680, is such that wh en a < 



25 



26 



27 



Qm, no 



29j for 



bound states can form in the potential. We refer the reader to Refs. 
detailed studies on the critical strengths of attractive central potentials, and for techniques 
to compute them. 

When the free energy is taken as potential term, the dissociation temperature of the 
[(n + l),£] state, denoted as T n i, is thus given by the solution of equation 

8/i 



m D {T nt 



3c/ 



nt ) 



(16) 



which is simply a rewriting of the definition (fl~5l) . It is worth mentioning that this result is 
independent of the explicit form of the strong coupling constant and of the screening mass. 
The only input, through the values of g n i, is the Yukawa form ([5]). Consequently, for another 
form of the radial potential like V(r,T) = —aa s (T) mfl(T) v[rriD{T)r}, with v(x) arbitrary, 
this relation would still be correct [a = 4/3 in Eq. fl!6p ]; only the value of g n £ would change. 
If we inject the definition (Q in (fIBl) . we obtain 

8/i 



T, 



(17) 



37 gni 

Since the model we developed here is physically relevant only above T c , and since g decreases 
with increasing T [see Eq. ffT5]) ]. bound states can form only if 



%lxa s (T) 



3m D (T) 



T=T C 



37 T c 



> goo- 



(18) 



It implies a lower bound on /x, which reads 



y, > -goolT c = 1.560 ± 0.036 GeV. 
8 



(19) 



Following the values of the Particle Data Book concerning the current quark masses 30[ , this 
lower bound only allows for bb bound states - for completeness, we recall indeed that the most 
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recent experimental values are m c = 1.25 ±0.09 GeV and mj = 4.45 ±0.32 GeV [30j. Let us 
note that, in a Schrddinger-based formalism, constituent quark masses are more commonly 
used. The intuitive idea is that a confined quark acquires an additional mass in a hadron 
because of the particles-antiparticles pairs which are created around it. The difference is 
particularly important for light quarks: Although the current mass of the u or d quarks is 



nearly zero, their constituent mass is often taken to be around 0.3 GeV 3JJ. The use of 
the constituent mass then allows to deal with a well-defined kinetic operator p 2 /2m. For 
heavy quarks, the current and constituent masses are essentially the same, since the current 
mass is much greater than the additional mass coming from any particles-antiparticles cloud. 
Nevertheless, the constraint (1191) forbids the existence of mesons made of at least one light 
quark above T c , either the current or the constituent quark mass is used. 

Following Eq. ffTTj) . the highest dissociation temperature is Too for a bb pair, thus the 
dissociation temperature of the IS bottomonium. It is readily computed that, in this case, 
Too — (1.426±0.070) x T c . The other values for T n £ are always smaller than T c even for the bb 
states, thus outside the validity range of our model. Interestingly, the prediction that even 
the heaviest mesons should be dissociated around 1.5 X T c is in qualitative agreement with a 
recent work based on the study of quarkonium correlators and spectral functions at nonzero 
temperature 32|. Let us note that, since our model does not take the spin interactions into 
account, we only have a single dissociation temperature for the both the rjb(lS) and T(1S) 
mesons. But, as it can be checked in Refs. jsl, Q, spin-dependent interactions are expected 
to affect the dissociation temperatures by less than 10 %. 



IV. INTERNAL ENERGY AS POTENTIAL TERM 

We turn now our attention to the case where the internal energy is taken as potential 
term in the eigenvalue problem to solve. Since, following Eq. ([9]), the screening mass is 
proportional to Ta s (T), the internal energy (T4J) can be rewritten as 

4 p-m D {T)r 

V(r,T) = U 1 (r,T) = -- [T a s (T)]' [1 - m D {T)r\ . (20) 

O T 

The sign of the internal energy crucially on the one of [T a s (T)]'. It is easily checked from 
definition ([5]) that, for T > T c , [T a s (T)] > 0. As an illustration, we plotted U^r, T) for 
T = 1.5 x T c in Fig. [3 
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0.02 




r (GeV) 1 



FIG. 2: Plot of Ui(r, T) computed for T = 1.5 x T c and with the optimal values of the parameters. 

Let us begin by a study of the possible bound states. The main equation (JT]) can again 
be rewritten in the form of Eq. (fT3"l) . but with 



P(x) 



x) 



and g > defined as 
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Sfi[Ta a (T)]' _ 8/i 
3 m D {T) 37 



x 
1 

T 



« S (T) 



(21) 



(22) 



The critical stren gth q nn is equal to 4.937 for the potential f[2~Tj) . We still refer the reader 



to Refs. 
exist if 



25 



26 



27 



28 



8fi 
37" 



29j for methods to compute g 00 . Consequently, bound states will 
8/x 



1 oATY 

T a s (T) 



37 T 



1 - 



ln(T/A CT ) 



> fi'oo- 



(23) 



Since the function of T in this last relation is ever decreasing for T > T c , a lower bound on 



fi is obtained when Eq. ( J23l) is computed for T = T c . One finds 



^ 37 T c 



hx/3 



8.213 ± 1.215 GeV. 



(24) 



_lnp + 1_ 

Apart from hypothetical bound states involving top quarks, this condition states that no 
bound state can survive after deconfinement provided the internal energy is used as potential 
term. The minimal reduced mass is indeed higher than any physically allowed reduced mass. 

If bound states cannot form at high temperatures, it should be however possible to find 
quasi-bound states in the potential Ui(r, T), i.e. resonance states with a positive eigenenergy 
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and a finite lifetime. This is due to the presence of a positive barrier in Ui(r, T), as it can be 
observed in Fig. [2j The function P(x) given by relation (]2ip has a minimum in x = M, with 
M = (1 + v^/2. The maximum value of the potential —gP(x) is thus reached in x = M, 
and quasi-bound states are expected to appear for 

£ < -gP(M). (25) 

These states will have a lifetime which goes smaller as £ gets closer of —gP(M). We stress 
that the results we discuss here concerning quasi-bound states are only valid in the s-channel. 
A rewriting of inequality (|25|) leads to an upper bound on the physical binding energy, that 
is 

4 e~ M 

e < 5 max (T) = -(M-l) — m D (T) [a s (T)T]\ (26) 

Figure [3] shows the evolution of the upper bound e max (T) with the temperature. Since it is 
always increasing, the formation of quasi-bound states seems favored at high temperatures. 
However, one can expect that, for increasing T, the average energy e of a qq pair in the 
deconfined medium will increase. Consequently, even if the particular shape of the internal 
energy gives more chance for the existence of quasi-bound states at high temperature, an 
ever decreasing number of qq pairs will have the appropriate energy e. A detailed study 
of these quasi-bound states thus requires a careful analysis involving statistical mechanical 
arguments, which is out of the scope of this paper. Finally, we can notice that the radius of 
these states should necessarily be smaller than x = M. It implies that their physical radius 
is such that r < M/m_o(T). This upper bound tends to zero as T becomes larger. 



V. SUMMARY OF THE RESULTS 



In this work, we have investigated the possible existence of mesons beyond the critical 
temperature. We have built a potential model from a commonly accepted Yukawa form for 
the free energy, that is Eq. ([2]). We used the well-known one- loop expression for the running 
coupling constant which is simpler than the two-loops expression and contains mainly the 



same physical _in 
result of Ref. 



ormation. For the screening mass however, we took the nonperturbative 



211 ]. which predicts that m D (T) oc a s (T)T instead of the usual dependence 



in a/« s (T) T. We have checked that the expression (jUJ) that we used for mc(T) is indeed 
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FIG. 3: Plot of the maximal energy e ma . x (T) for the existence of quasi-bound states states as a 
function of the temperature in units of T c . 



in better agreement with recent lattice QCD simulations ll| than the perturbative-like 
form (fTUl) in the temperature range (1 — 3) x T c . 

We provide in this paper a simple method, which allows us to get mainly analytical 
results and can be applied to different versions of the model discussed here. In particular, 
we have shown that the determination of the dissociation temperature is ruled by a single 
dimensionless parameter which corresponds to the strength of the potential in a rescaled set 
of coordinates. Constrains on the dissociation temperature are thus obtained from constrains 
on the strength of the potential such as it is attractive enough to admit bound states. We 
observed very different behaviors following that the free energy or the internal energy is 
chosen as potential term in the Schrodinger equation we study. 

When the free energy is used, we have shown that the dissociation temperature of a 
given meson, that is the temperature at which it becomes unbound, can be analytically 
obtained. Logically, the dissociation temperature decreases for increasing n and £, but also 
for decreasing \x. We deduced from Eq. (TPTj) that the only state which could be present 
above T c is the IS bottomonium, with a dissociation temperature of (1.426 ± 0.070) x T c , 



in qualitative agreement with another recent work 32]. 

Using the internal energy as potential term leads to a rather different picture. The 
existence of bound states is still only ruled by a dimensionless parameter [see Eq. (1221) ]. but 
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we have shown that "true mesons", i.e. qq bound states, cannot exist anymore above T c 
in this case. However, it is important to notice that the potential now exhibits a positive 
barrier even at zero angular momentum. This allows for the existence of quasi-bound states 
(metastable states) in the continuum, i.e. for positive eigenenergy. In this picture, a few 
heavy qq pairs could thus subsist in the medium above T c as compact (small radius) quasi- 
bound states. 
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